AIAA-2000-1660 


BENDING BOUNDARY LAYERS IN LAMINATED-COMPOSITE 
CIRCULAR CYLINDRICAL SHELLS* 

Michael P. Nemeth* and Stanley S. Smeltzer, III** 

Mechanics and Durability Branch, NASA Langley Research Center 
Hampton, Virginia 2368 1-2199 


A b str ag j 

An analytical, parametric study of the attenuation of 
bending boundary layers or edge effects in balanced and 
unbalanced, symmetrically and unsymmetrically laminat- 
ed thin cylindrical shells is presented for nine contempo- 
rary material systems. The analysis is based on the linear 
Sanders-Koiter shell equations and specializations to the 
Love-Kirchhoff shell equations and Donnell’s equations 
are included. Two nondimensional parameters are inden- 
tified that characterize and quantify the effects of laminate 
orthotropy and laminate anisotropy on the bending bound- 
ary-layer decay length in a very general and encompassing 
manner. 

A substantial number of structural design technology 
results are presented for a wide range of laminated-com- 
posite cylinders. For all the laminate constructions con- 
sidered, the results show that the differences between 
results that were obtained with the Sanders-Koiter shell 
equations, the Love-Kirchhoff shell equations, and Don- 
nell’s equations are negligible. The results also show that 
the effect of anisotropy in the form of coupling between 
pure bending and twisting has a neglible effect on the size 
of the bending boundary-layer decay length of the bal- 
anced, symmetrically laminated cylinders considered. 
Moreover, the results show that coupling between the var- 
ious types of shell anisotropies has a negligible effect on 
the calculation of the bending boundary-layer decay 
length in most cases. The results also show that in some 
cases neglecting the shell anisotropy results in underesti- 
mating the bending boundary-layer decay length and in 
other cases it results in an overestimation. 
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laminate membrane stiffnesses 


modifed laminate stiffnesses defined in 
Appendix B 

laminate membrane-bending coupling 
stiffnesses 

modifed laminate stiffnesses defined in 
Appendix B 

first-order correction factor for 
anisotropy parameter 
attenuation or decay lengths 
laminate bending stiffnesses 


stiffness coefficients 
major and minor principal lamina 
moduli, respectively 
lamina shear modulus 
cylinder wall thickness and length 
axial and circumferential bending and 
twisting stress resultants, respectively 
axial, circumferential, and shear 
membrane stress resultants, respectively 
nondimensional orthotropy parameter 
loading function appearing in bending 
boundary-layer differential equation 
transverse-shear stress resultants 
cylinder radius 

constant coefficients of bending 
boundary-layer differential equation 

modified shear stress resultant defined in 
Appendix B 
strain-energy density 
normal-displacement component 
axial coordinate of cylinder 
attenuation-length tolerance parameter 

axial, circumferential, and shear 
membrane strains, respectively 
circumferential, angular coordinate 

axial and circumferential bending and 
twisting strains, respectively 
constant defining different shell theories 
lamina major Poisson’s ratio 
lamina fiber angle 
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Introduction 

The term, "bending boundary layer," refers to lo- 
calized zones of bending stresses and deformations that 
appear in practically every type of thin shell structure. 1 
Bending boundary layers are caused by edge support 
conditions; by localized mechanical loads, heating, or 
cooling; and by abrupt changes in stiffness, such as that 
caused by a cutout, a crack, or a stiffener. All of these 
effects may be real concerns in a given preliminary de- 
sign for an advanced aerospace vehicle made of laminat- 
ed-composite materials. Thus, it is useful to have 
nondimensional parameters that characterize the effects 
of shell geometry and laminate orthotropy and anisotro- 
py on the extent of bending boundary layers and that can 
be used to help guide the development of a design. For 
example, an optimal design for a pressure vessel might 
be one that exploits the membrane load-carrying action 
of a shell and minimizes zones of local bending stresses. 
Meaningful nondimensional parameters could be used at 
the preliminary design stage to identify families of lami- 
nates and material systems that exhibit relatively small 
bending boundary layers. Moreover, a meaningful esti- 
mate of the size of a bending boundary layer in a shell is 
very useful for determining an adequate First-approxima- 
tion Finite-element model for a complex shell structure. 
Without a proper understanding of the extent of a bend- 
ing boundary layer, it is possible to have a Finite-element 
model that could miss a significant part of the structural 
response in a region where failures are often initiated by 
high interlaminar stresses. Furthermore, apriori knowl- 
edge of the extent of bending boundary layers is useful in 
determining the instrumentation locations in structural 
verification tests or in material characterization tests. 2 In 
addition, knowledge of how laminate construction af- 
fects the extent of a bending boundary layer is useful for 
understanding how nonlinear prebuckling deformations 
affect the buckling behavior of cylindrical shells. 

Studies of the behavior of axisymmetric, bending 
boundary layers in right-circular, cylindrical shell struc- 
tures made of orthotropic or anisotropic materials and 
with finite length have been presented, to at least some 
extent, in Refs. I through 13. In the discussion that fol- 
lows, reference is made to unbalanced and balanced lam- 
inates that are either symmetrically or unsymetrically 
laminated. Herein, the term unbalanced laminate is used 
to indicate that coupling between pure extension or con- 
traction and shearing is present in a laminate. The term 
unsymmetric laminate is used to indicate coupling be- 
tween any of the components of bending action with any 
of the components of membrane action. A fully anisotro- 
pic laminate would include both of these types of anisot- 
ropy in addition to the anisotropy that is manifested by 
coupling between pure bending and twisting action that 


is sometimes exhibited by balanced, symmetric lami- 
nates. 

In Ref. 1 , an analysis is presented and an expression 
for the attenuation or decay length of the bending bound- 
ary layer for a specially orthotropic cylinder that is sub- 
jected to edges loads, internal pressure, and heating is 
given. These equations, and the accompanying results, 
are based on the linear Love-Kirchhoff shell equations. 
In Ref. 3, an analytical solution that is based on Don- 
nell’s simplifications to the linear Love-Kirchhoff shell 
equations is given for fully anisotropic cylinders that are 
subjected to lateral pressure and edge loads. Results that 
show the effect of laminate anisotropy on the edge mo- 
ment are presented for a clamped two-ply shell that is 
subjected to internal pressure. In addition, a discussion 
is presented that suggests that solutions that are based on 
Donnell’s equations should be accurate for laminates 
that are not highly anisotropic. An analytical study of 
bending boundary layers in unbalanced, symmetrically 
laminated cylinders, that is also based on Donnell’s 
equations, is presented in Ref. 2. The aim of this study 
was to determine a suitable gage section in a laminated- 
composite tube that is to be used for a material character- 
ization test. Results are presented for unidirectional, he- 
lical-wound tubes. 

An analytical solution for bending boundary layers 
in unbalanced, symmetrically laminated and balanced, 
unsymmetrically laminated circular cylindrical shells 
that are subjected to internal pressure and thermal loads 
is presented in Ref. 4. The solution is also based on Don- 
nell’s linear equations and numerical results are present- 
ed for filament-wound cylinders made of heat-treated 
carbon-carbon material. A study that focuses mainly on 
prebuckling deformations, with bending boundary lay- 
ers, in homogeneous, orthotropic and unsymmetrically 
laminated cross-ply cylinders that are subjected to axial- 
compression loads and lateral pressure loads is presented 
in Ref. 5. The effects of the bending boundary layers on 
the buckling response are examined for several laminate 
constructions, but the general effects of the laminate con- 
struction on the extent the boundary layers are not dis- 
cussed. 

A pair of complex conjugate, fourth-order equa- 
tions that are based on Flugge’s corresponding 
equations 14 and that can be solved in closed form are de- 
rived for specially orthotropic, circular cylindrical shells 
in Ref. 6. Moreover, eigenfunction solutions are present- 
ed that include the solution for the axisymmetric bending 
boundary layer and several simplifed equations are pre- 
sented and their relative accuracy is analyzed. In Ref. 7, 
a study of bending boundary layers in transversely iso- 
tropic circular cylindrical shells is presented. This study 
examines the attenuation characteristics of bending 
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boundary layers by applying an asymptotic method to the 
linear, three-dimensional elasticity equations, and pre- 
sents ordcr-of-magnitude estimates for the stresses and 
displacements for a wide range of ratios of the two prin- 
cipal elastic moduli, [n Ref. 8, an analytical solution for 
an unbalanced, unsymmetrically laminated circular cy- 
lindrical shell that is subjected to internal pressure is pre- 
sented that is based on a variant of the Love-Kirchhoff 
shell theory, which uses an expression for the change in 
surface twist that was given by Timoshenko. Numerical 
results are also presented for a two-ply shell that demon- 
strate the coupling effects of the shell anisotropies. 

The bending boundary layers of an unbalanced, un- 
symmetrically laminated circular cylindrical shell that is 
subjected to axial compression, torsion, or thermal load- 
ing are investigated in Ref. 9. Results are also presented 
that demonstrate the coupling effects of the shell 
anisotropies. In addition, results are presented for two 
more conventional unsymmetric laminates and a typical 
quasi-isotropic laminate. In Refs. 10 and 11, bending 
boundary layers are also examined for balanced, sym- 
metrically laminated and balanced, unsymmetrically 
laminated cylindrical shells, in the context of nonlinear 
prebuckling deformations that occur as a result of com- 
pression and thermal loads. In particular, the effects of 
laminate stacking sequence on the extent and character 
of the bending boundary layers are presented for two 
groups of three similar laminates. Two of the laminates 
are unsymmetric. In Ref. 1 2, a linear analysis is present- 
ed that focuses mainly on balanced, symmetrically lami- 
nated cylinders, and an expression is given for the length 
of the bending boundary layers near the cylinder ends 
that is based on the Love-Kirchhoff shell equations. 

Most recently, Golden veizer’s static-geometric du- 
ality principle 15 has been used in Ref. 13 to reduce the 
Sanders-Koiter equations 16,17 for fully anisotropic, right- 
circular cylindrical shells to two coupled fourth-order 
equations that use a stress and a curvature function as the 
unknown, primary field variables. The reduction is done 
by adding certain negligibly small terms to the stress- 
strain relations, which are intrinsically in error because 
they must be established experimentally. The approach 
demonstrates how the static-geometric duality principle 
can be used to reduce greatly the amount of algebra need- 
ed to obtain results. Eigenfunction solutions are also pre- 
sented for specially orthotropic cylinders that are in 
agreement with corresponding results presented in Ref. 
6. Moreover, asymptotic formulas that can be used to de- 
termine axisymmetric bending boundary layer attenua- 
tion lengths and the decay of other unsymmetric, self- 
equilibrated edges loads are given. 

With the exception of Ref. 13, explicit expressions 
for estimating the size of axisymmetric bending bound- 


ary layers in fully anisotropic, right-circular cylinders arc 
not found in the literature. Moreover, there appears to be 
even fewer results for laminated-composite shells made 
of contemporary material systems and essentially no 
substantial parametric studies. The present paper focus- 
es on developing meaningful estimates of attenuation 
lengths of bending boundary layers in balanced and un- 
balanced, symmetrically and unsymmetrically laminated 
circular cylinders. The analysis is based on the linear 
Sanders-Koiter shell equations and contains the Love- 
Kirchhoff shell equations 1 and Donnell’s equations 1 as 
special cases, and is somewhat similar to the analyses 
presented by Reuter 3 4 and Chaudhuri, et. al. 8 With these 
equations, explicit expressions are obtained and nondi- 
mensional parameters are presented that characterize the 
effects of cylinder geometry and laminate construction 
on the size of a bending boundary layer in a very general 
manner. In particular, generic design curves are present- 
ed that use the nondimensional parameters to show the 
effects of laminate orthotropy and anisotropy on the at- 
tenuation length in a concise and encompassing manner. 
In addition, values of these parameters are presented for 
a very wide range of orthotropic and anisotropic laminate 
constructions. Also, differences in the results that were 
obtained in the present study by using the Sanders-Koiter 
shell equations, the Love-Kirchhoff shell equations, and 
Donnell’s equations are discussed. 

Analysis 

The ordinary differential equation that governs the 
axisymmetric bending behavior of a right-circular cylin- 
der that is subjected to edge loads or displacements and 
surface tractions is obtained by first specializing the lin- 
ear Sanders-Koiter shell equations, that are given in Ap- 
pendix A, for axial symmetry. For the equations 
presented herein, x and 0 denote the axial and circum- 
ferential coordinates of a right-circular cylinder, respec- 
tively, and the specialization to axial symmetry is 
obtained by eliminating all terms in the equations that are 
differentiatied with respect to the circumferential coordi- 
nate, 0. The resulting set of equations for axisymmetric 
behavior are given in Appendix B. The ordinary differ- 
ential equation that governs the axisymmetric bending 
behavior of a right-circular cylinder that is subjected to 
edge loads or displacements and surface tractions is de- 
rived in Appendix B and is given by 

dF +4S & +4Qw=P(x) (1) 

where S, Q, and P(x) are defined in Appendix B by Eqs. 
(B55), (B56), and (B57), respectively, and w(x) is the 
radial deflection that is positive-valued when outward. 
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The coefficients of Eq. (1) depend on the subscripted A, 
B, and D constitutive terms of classical Love-Kirchhoff- 
type laminated shell theory (e.g., see Ref. 18, pp. 1 90- 
202) and the radius of the cylinder middle surface, R. 

To determine the specific form of the solution to 
Eq. (1), it is useful to examine the positive-definiteness 
conditions on the strain-energy density function. The 
strain energy density function for this problem is given 
by 


constitutive matrix in Eq. (4) be positive-valued. More- 
over, by rearranging the strain energy density function 
into the form 




A,, A l6 A l2 B lt 
A |6 S 66 A :6 b 16 
A 12 A 26 A 22 ® 12 

B ii B 16 B }2 D , , 



(7) 


rU = N X + N 6 £ 0 + N xe ^ 0 + MX + M X0 K X0 (2) 


the following additional positive-definiteness condition 
is obtained 


where N x , N 0 , and N x0 are the membrane stress result- 
ants; M x and M xe are bending stress resultants; £*, Eq, 
and y° x9 are the middle-surface membrane strains; and 

K° and K x0 are are middle-surface bending strains. By 
using Eqs. (B22), (B23), and (B28), the strain-energy 
density function is expressed as 


A,, A l6 
A 16 ^66 


“ A 1^66 ~ A 16 > 0 


( 8 ) 


The homogeneous solution for Eq. (1) involves the 
square root of the quantity Q - S 2 . By using Eqs. (B55) 
and (B56), this quantity is given by 


2 11 = NX + NX + Tyi + MX (3) 


Q-S 2 = 


4C x C 3 -C 2 
\6C] 


(9) 


The strain energy density is expressed in terms of the 
strains and constitutive terms by using the constitutive 
equation given by Eq. (B29); that is, 


Substituting Eqs. (B41) - (B43) into Eq. (9) and simpli 
fying, the quantity Q - S 2 is found to be given by 




1 

2 



A u A 12 S 16 B n 
A l2 A 22 ^26 Bj 2 
A 16 A 26 A 6 6 ® 16 

B n B i2 B l6 Du 



(4) 


Q-S 


2 __ 1 

4C; 


An A 12 A j 6 B u 
A 12 A 22 A 26 B 12 
A 16 A 26 A ^ B l6 
B n B 12 B 16 D u 


( 10 ) 


The stiffness terms in Eq. (4) that have overbars are 
defined by Eqs. (B31) - (B35) and are functions of the 
shell wall thickness-to-radius parameter, h/R. By 
enforcing positive definiteness of the strain energy den- 
sity function (e.g., see Ref. 19), the requirement that the 

diagonal terms A n> A 22 , A^, and D n be positive-val- 
ued is obtained. Moreover, the following determinants 
are positive valued 

= A ijA 2 2 — A 12 > 0 (5) 


A u A 12 
A 12 A 22 


An A i2 A l6 
A ,2 A 22 A 26 
A 16 A 26 A 66 



22 



66 ~ 


It follows logically, that Q - S 2 > 0 because the posi- 
tive-definiteness of the strain energy density function 
requires that the determinant in Eq. (10) be positive val- 
ued. Moreover, Q - S 2 > 0 implies that Q > 0, and Q > 

Q 

0 implies that pr > 0 . Equations (6), (8), and (B43) 

Q 

indicate that C 3 > 0. Thus, pr > 0 yields the condition 

c i 

that C { = D n e> 0 (see Eqs. (B49) and (B50)). Because 
D n > 0, e > 0. To enunciate the positive valuedness of 
Q, it is convenient to introduce the expression 


T 2 = Q = 


1 


4R‘a 22 D n e 


( 11 ) 


such that T 2 - S 2 > 0, and to express Eq. (1) as 


A i|A 26 - A 22 A 2 6 + 2A 12 A i 6 A 26 > 0 (6) 

Likewise, positive definiteness of the strain energy den- 
sity function also requires that the determinant of the 


d^ +4S d^ +4T:w=P(x) < 12) 


Equation (12) is a linear, fourth-order, nonhomoge- 
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ncous ordinary differential equation with constant coef- 
ficients. The characteristic equation of Eq. (12) is given 
by 

V + 4S X 2 + 4T 2 = 0 (13) 

Using the knowledge that T 2 - S 2 > 0, the roots of the 
characteristic equation are obtained from the quadratic 
formula; that is, 

(^) 12 = 2(-S±«Vt 2 -S 2 ) (14) 

where i - . Solution of this equation for X yields 

four roots of Eq. (13) that are pairs of complex conju- 
gates that are given by 

^i.2,3,4 ~ ~ ^ ^ ^ (15) 


Formulas for the Attenuation Len gth 

Formulas for the attenuation or decay length of the 
bending boundary layers arc obtained by first noting that 
the response quantities for the region near x = 0 are 

bounded by the two functions i F|C ^ T_S * and that the 
response quantities for the region near x = L are bounded 

by the two functions i F 3 e ” s (L Let d denote the 
length for which the solution attenuates or decays to a 
value of e times the amplitude F t or F 3 . A reasonable es- 
timate of the attenuation length or decay length d is ob- 
tained by replacing x and (L - x) with d in the 
exponential terms of Eq. (16b), and by noting that the 
amplitude of w(x) is attenuated by the exponential terms. 
Thus, the attenuation length or decay length is given by 

e -*/ T -s</ = e w hi c h yields 

d = -Ine (T - S)~ l/2 < j (17) 


The homogeneous solution of Eq. (12) can be writ- 
ten as follows 

w H (x) = sin[yT + Sx + K 2 ] + 


By using Eqs. (B55) and (1 1), Eq. (17) is expressed as 


d = d° 

vT*h fRh 


(18) 


K 3 e-^<L-’> s in[,/TTTx + K 4 ] (i 6a) 

where x e [0,L] . The symbols K,, K 2 , K 3 , and K 4 are 
real-valued constants that are determined from the 
boundary conditions given by Eqs. (B18) and (B19). 
The solution given by Eq. (16a) represents a damped, 
oscillatory response that decays from each end of the 
cylinder. The regions near the edges of the cylinder, 
where the amplitude of w H (x) is the largest are called the 
bending boundary layers. All response quantities that 
exhibit bending boundary layers involve derivatives of 
Eq. (16a) and can be expressed in the general form 

F(x) = F l e” v ^ r s* sin[v/T + S x + F 2 ] + 

F 3 e _ '^ a '‘ , sin[vTTs x + F 4 ] (16b) 

where F, through F 4 are constants. 

When the length of the bending boundary layers are 
less than half of the cylinder length, which is typical, 
Eqs. (16) can be partitioned into one part that applies to 
the edge x = 0 and the other that applies to the edge x = 
L. The response quantities for the region near x = 0 are 
obtained by setting F 3 = 0 in Eqs. (16). Similarly, the 
response quantities for the region near x = L are obtained 
by setting F, = 0 in Eqs. (16). 


where d° is the attenuation length, in which anisot- 
ropy is neglected, that is given in nondimensional form 
by 


d° 

VRh 



(19) 


The symbol h is the shell wall thickness, and 0 and /! 
are nondimensional orthotropy and anisotropy parame- 
ters or factors, respectively, that are given by 


0 = 


12 A ||Dn 
(A n A 22 -A 2 2 )h 2 


(20a) 




(A n A 22 - A 12 ) 
An 


<*22^ 


\j a 2 2 D i 


(20b) 


i c 


where the symbols in these equations are defined in 
Appendix B. 

Other useful forms of Eq. (20a) are obtained by in- 
troducing an effective membrane Poisson’s ratio 
A 

v m - i which is the geometric mean of the two 

V A 1 1 A 22 

Poisson effects associated with the inplane principal di- 
rection of a homogenized orthotropic material. By using 
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this effective membrane Poisson’s ratio, Eq. (20a) is ex- 
pressed as 


the only anisotropic constitutive terms are D I(S and D 2fi , 
and the anisotropy factor is given by A = Ve where 


0 = 

12 D u 

1/4 

(21) 


( h V 

A 22 h 2 (l -v 2 ) 


(i) 

For a single-layer of homogeneous, specially orthotro- 

A«,D u h 2 

[ p Waj.’ 


(25) 


pic material, v m = V v i2 v 2i » = 1 » and 

E, 


0 = 


E 2 (l-v l2 v 21 ) 


( 22 ) 


which, when substituted into Eq. (19), yields results 
identical to the results presented by Kraus 1 , where the 
decay tolerance is given by e = e‘\ Likewise, for a sin- 
gle-layer of isotropic material with an arbitrary thick- 
ness, v m = v, A = 1, and 


For thin-shell theories, such as the Sanders-Koiter the- 
ory and the Love-Kirchhoff theory, -20- This result 

suggests that a useful approximation to Eq. (25) and the 
anisotropy factor can be obtained from a power series 

h 

expansion for small values of This process yields 




D;, 


4 1RJ AJ) u h 2 


(26) 


0 = 


1 

Vd-v 2 ) 


(23) 


A 90%-decay length (e = .1) that is a good approxima- 
tion to the behavior of homogeneous, metallic shells is 
d° f~h~ 

given by ^ - 1-79 y ^ . Applying this formula, for 

example, to the Space Shuttle solid rocket booster 
described in Refs. 20 or 21 (R = 72 in., h= 0.5 in.) gives 

d° =0.15 R= 10.8 in. 

It is interesting to note that the differences between 
the attenuation lengths that are based on the Sanders- 
Koiter, the Love-Kirchhoff, and Donnell’s equations ap- 
pear in the coefficient e and in the symbols with over- 
bars in Eq. (20b) for the anisotropy factor A (see Eqs. 
(B22), (B3 1) - (B35), and (B45) - (B49). For these equa- 

3 

tions, the Sanders-Koiter theory is given by P - j and 

the Love-Kirchhoff theory is given by p = 1. Donnell’s 
equations are given by p = 0. For isotropic and specially 
orthotropic cylinders, A - 1 and the three sets of shell 
equations yield identical results. Similarly, for antisym- 
metric cross-ply cylinders (A 16 = A 26 = D I6 = D 26 = B l6 = 
B 26 = B l2 = B 66 = 0) 



and the three shell theories yield identical results. 


D 


In this expression, 0 < p < ^ and 0 < — — - 


2 < 1 . 


Thus, the approximate formula for A indicates that for 
most practical applications of thin-shell theory, the dif- 
ferences between the three different shell theories con- 
sidered herein, and the effect of the flexural anisotropy 
of a general symmetrically laminated cylinder, are negli- 
gible. 

A simplified formula for the anisotropy factor can 
be derived for the general expression for A that is given 
by Eq. (20b). For this case, the following power series 

expansions for small values of g are used 


e=e 0 + ei(^) + e 2 (^)'+ ■ (27) 

&i2 = a i2 + ^ u2 (r) + ^ 2,2 (r) (28) 

^22 = a 22 + ^>22^ p[ j + a 222^J^ j + ” (29) 

^26 = (30) 

^21 = b 2 | + b I2t|^ j + + *' (31) 


Substituting Eqs. (27) - (31) into Eq. (20b) and expand- 
ing the resulting expression in a similar manner yields 


Simplified Formulas for A 

For balanced, symmetrically laminated cylinders, 


* ~* 0 + *{r) + *{r) + ' 


(32) 
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The coefficient /i l is a very complicated expression, and 
as a result, the following first-order approximation of/^ 
is used herein; that is, 





(33) 


where /^ 0 is the value of Eq. (20b) with ]X ~ 0, which is 
the anisotropy factor that corresponds to the use of Don- 
nell’s equations. This expression is given by 


^0 = 


(A^h-A’,) . 

1/4 

1 bii 

A d 22^0 

n II 


^ a 12 D„ c c 


(34) 


In this expression, a^ and b 2l are obtained from Eqs. 
(B45) - (B48) by setting pi = 0 in Eqs. (B31) - (B35). 
The expression for e 0 is obtained from Eq. (B49) in a 
similar manner. The term represents a first-order 
correction to the results that correspond to Donnell’s 
equations and is given by 


\J ^ 22 ^ 11^0 | 

a^, + a t22 tf 0 

| + 2a :: j 

b , 2|£ 0 1 

bj*^i 

| 2a, 22 b 2J e 0 

4a 2 * 0 

\j a : A 

1*0 ^ 21 



(35) 


where the terms that appear in Eq. (35) are given in 
Appendix C. In addition, further simplifications to 
and C\ are also presented in Appendix C for unbalanced 
and balanced, symmetric laminates and for balanced, 
unsymmetric laminates, that include the subclasses of 
general antisymmetric laminates, antisymmetric cross- 
ply laminates, and antisymmetric angle-ply laminates. 
The relative size of and its contribution to Eq. (33) 
are examined parametrically in the subsequent section 
of the present study. 


Results and Discussion 


Equations (18) and (19) form the basis for the para- 
metric study presented herein. In particular, the two 
equations isolate the contributions of shell orthotropy 
and shell anisotropy to the bending boundary-layer de- 
cay length with nondimensional parameters and imply 
the generic design-chart representations that are illustrat- 
ed in Figs. 1 and 2. In Fig. 1 , generic results are present- 
ed that show the nondimensional, 90%-decay length 


given by 


</Rh 


e-o.i 


as a function of the orthotropy pa- 


rameter 0 , for selected values of the anisotropy parame- 
ter/^. A 90%-decay length was selected herein to yield 
an accuracy that is approximately to within the accuracy 
of the experimentally determined material properties, but 


other values could be used. 

In a manner similar to Fig. I , Fig. 2 shows the non- 
dimensional, 90%-decay length as a function of the 
anisotropy parameter/^, for selected values of the orthot- 
ropy parameter 0 . Results that correspond to balanced, 
symmetrically laminated cylinders are given by a value 
of A = l and results that correspond to an isotropic shell 
wall are indicated in the figures by the filled circle with 
an ordinate value of 1 .79. Overall, these two figures rep- 
resent results that are applicable to a vast range of lami- 
nate constructions, and provide a common basis for 
comparison of regular and hybrid laminates made of dif- 
ferent material systems and laminate stacking sequences. 
In general, the figures show increases in the nondimen- 
sional 90%-decay length with increases in either of the 
orthotropy parameter 0 or the anisotropy parameter /t. 
In addition, the results in Figs. 1 and 2 clearly indicate 
the effect of neglecting shell-wall anisotropy on the at- 
tenuation length of a bending boundary layer. 

The actual value of the nondimensional, 90%-de- 
cay length depends on the particular values of the orthot- 
ropy and anisotropy parameters of a given laminate. 
Thus, additional results are presented subsequently that 
show how the orthotropy parameter 0 and the anisotropy 
parameter/^ vary with laminate construction. In particu- 
lar, values of 0 and are presented first for balanced 
and unbalanced symmetrically laminated cylinders. 
Then, values are presented for balanced and unbalanced 
unsymmetrically laminated cylinders. Nine different 
contemporary material systems were used to generate 
these results. These material systems include boron-alu- 
minum, S-glass-epoxy, a typical boron-epoxy, AS4/ 
3501-6 graphite-epoxy, AS4/3502 graphite-epoxy, IM7/ 
5260 graphite-bismaleimide, Kevlar 49-epoxy, IM7/ 
PETI-5, and P-100/3502 pitch-epoxy materials. The me- 
chanical properties of these material systems are present- 
ed in Table 1 and the nominal ply thickness that was used 
is 0.005 in. 


Balanced. Symmetrically Laminated Cylinders 

Symmetrically laminated shell walls are character- 
ized mathematically by values of zero for the subscripted 
B terms that appear in the constitutive equation, Eq. 
(A 15). In addition, balanced, symmetrically laminated 
shell walls do not exhibit coupling between extension 
and shear, which is characterized by A sfi = A 26 = 0 in Eq. 
(A 15). Shell walls of this class are strictly specially 
orthotropic for many laminates. However, for some wall 
constructions, balanced, symmetric laminates exhibit 
anisotropy in the form of coupling between pure bending 
and twisting of the shell wall. This type of anisotropy is 
manifested by nonzero values of the D i6 and D :6 consti- 
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tutivc terms in Eq. (A 15). However, the discussion of 
Eq. (26) that has been given herein indicates that this 
type of anisotropy is negligible for thin shells and that the 
differences between results obtained from the Sanders- 
Koiter, the Love-Kirchhoff, and Donnell theories are in- 
significant. Moreover, 1 for this class of laminated- 
composite shell walls, and the attenuation behavior is 
governed by the nondimensional orthotropy parameter 0 
that is given by Eq. (20a). Furthermore, Eqs. (18) and 
(19) indicate that the attenuation length is a constant 
multiple of the orthotropy parameter that depends on the 
attenuation-tolerance parameter £. For this case, trends 
that are exhibited by 0 are identical to those exhibited by 
the attenuation length based on any value of £. 

Values of the orthotropy parameter 0 are presented 
in Fig. 3 for single-ply, homogeneous, specially ortho- 
tropic and isotropic shell walls, with arbitrary thickness, 
as a function of the ratio of the principal elastic moduli, 
E 2 /E r For these results, the orthotropy parameter is giv- 
en by Eq. (22) and is expressed in the following more 
convenient form 



One curve, that is essentially several coincident curves, 
is shown in the figure that corresponds to general results 

for 0.2 < v l2 ^ 0.35 . In addition, specific results for the 
nine material systems considered herein and for a typi- 
cal aluminum and a steel are indicated in the figure by 
the square symbols. The results in Fig. 3 indicate that 
the effect of variations in the major Poisson’s ratio on 
the orthotropy parameter 0 are small compared to the 
effect of variations in the ratio of the principal elastic 
moduli. Moreover, the results show that 0 decreases 
rapidly as the ratio of the principal elastic moduli 
increases, particularly for values of E 2 /E t less than 
approximately 0.1, which corresponds to most of the 
contemporary orthotropic materials considered herein. 
Figure 3 also shows that an isotropic material corre- 
sponds to 0 ~ 1 . 

Values of the orthotropy parameter 0 for the sin- 
gle-ply, homogeneous, specially orthotropic cylinders 
investigated by Cheng and He 6 were also obtained. A 
comparison of the results obtained in the present study 
with the corresponding results of Ref. 6 are presented in 
Table 2 for boron-epoxy, glass-epoxy, and graphite-ep- 
oxy materials and for the cylinder radius-to-thickness ra- 
tio R/h = 208.3 1 1 . Moreover, a range of results is shown 
for Ref. 6 which corresponds to various simplifications 
that were used in the equations that govern the response. 


The actual material properties that were used arc given in 
Ref. 6. In this table, the quantity used for comparison is 
given by 



(37) 


which is the real part of the exponent p that appears in 
the eigenfunction solution used by Cheng and Ho (n = 0 
in Eq. (25) of Ref. 6; see also Eq. (47) of Ref. 1 3), which 
corresponds to the decay or attenuation of the response. 
The orthotropy parameter shown in Eq. (37) is defined 
by Eq. (36). The results in Table 2, show very good 
agreement (less than 1% difference) for all three materi- 
als. In addition, the results obtained herein that are 
shown in Table 2 for the boron-epoxy material are also 
in excellent agreement with the corresponding results 
presented by McDevitt and Simmonds. 13 

Values of the orthotropy parameter 0 are presented 
in Fig. 4 for multilayered [(±<J>) m ] s laminates made from 
the nine material systems as a function of the Fiber angle 
<)>, which is measured from the x-axis toward the 0-axis. 
The results are independent of the stacking sequence 
number m and show a wide variation in 6 with the ma- 
terial system. The results also show, for the most part, a 
wide variation in O with the Fiber angle <|> and show a re- 
duction in 0 as the Fiber angle increases from zero to 
ninety degrees. The largest value (2.93) and the smallest 
value (0.34) of O are exhibited by the unidirectional 
laminates made from P-100/3502 pitch-epoxy material, 


equal to 5.13 and 


e»o.t 


and correspond to values of 

0.59, respectively. Moreover, the greatest variation in O 
with the Fiber angle (approximately 8.7 times) is exhibit- 
ed by the laminates made from P-100/3502 pitch-epoxy 
material. The smallest variation is exhibited by the lam- 
inates made from the boron-aluminum material. 

Results are presented in Fig. 5 that show the values 
of the orthotropy parameter for [(±45/0 2 )J s , [(0 2 /±45) m ] s , 
[(±45/90 2 )J s , [(90 2 /±45)J s> [(±45/0/90 ) J s and [(0/90/ 
±45) m ] s laminates made of IM7/5260 graphite-bismale- 
imide material for values of the stacking sequence num- 
ber m = 1 to 6. Values of O range from approximately 
1.53 to 0.64. These results show that the curves for the 
[(±45/0 2 )J s and [(0 2 /±45)J s laminates approach O ~ 

1 .41 as m increases to a value of 6, with the curve for the 
[(0 2 /±45) m ] s laminates converging from above and the 
other curve converging from below. The higher values 
of 0 for the [(0 2 /±45) m ] s laminates, are attributed to the 
higher axial bending stiffness that is obtained by placing 
the zero-degree plies at the outer surfaces of the lami- 
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nates, particularly, for the lower values of the stacking 
sequence number m. Similarly, the results in Fig. 5 show 
that the curves for the [(±45/90 2 )J s and [(90 2 /±45) m ] s 
laminates approach 0 ~ 0.76 as m increases to a value 
of 6, with the curve for the [(±45/90 2 )J s laminates con- 
verging from above and the other curve converging from 
below. Likewise, the results in Fig. 5 show that the 
curves for the [(±45/0/90) m ] s and [(0/90/±45)J s quasi- 
isotropic laminates approach 0~ 1.03 as m increases 
to a value of 6, with the curve for the [ (0/90/ ±45) m ] s lam- 
inates converging from above and the other curve con- 
verging from below. 

Overall, the results in Fig. 5 indicate that the [(±45/ 
0 2 ) m ] s and [(0 2 /±45)J s laminates exhibit higher values of 
the orthotropy parameter than the [(±45/0/90) m ] s and [(0/ 
90/±45) m ] s quasi-isotropic laminates, which exhibit 
higher values of the orthotropy parameter than the [(±45/ 
90 2 )J S and [(90 2 /±45)J s laminates. This trend corre- 
sponds to a reduction in the value of 0 as the axial bend- 
ing and extensional stiffnesses of the laminates decrease. 

Results are presented in Fig. 6 that show the effect 
of the nine material systems considered herein on the 
orthotropy parameter for the [(0 2 /±45) m ] s laminates. 
Values of 0 range from approximately 1.67 for P-100/ 
3502 pitch-epoxy material to 1.09 for boron-aluminum 
material. Most of the materials exhibit values of 0 in the 
range of approximately 1 .4 to 1 .6. All of the curves show 
about the same reduction in 0 as the stacking sequence 
number m increases. 

Results similar to those in Fig. 6 are presented in 
Fig. 7 that show the effect of the nine material systems 
on the orthotropy parameter for the [(±45/0/90)J s and 
[(0/90/±45) m ] s quasi-isotropic laminates. These results 
show a much smaller variation in the orthotropy param- 
eter with material system and stacking sequence number 
for the quasi-isotropic laminates than for the [(0 2 /±45)J s 
laminates in Fig. 6. In particular, values of 0 for the qua- 
si-isotropic laminates range from approximately 1.15 to 
1 . The largest values of O in Fig. 7 are exhibited by lam- 
inates from P-100/3502 pitch-epoxy material. More- 
over, the results show a larger variation in 0 with 
stacking sequence number for the [(0/90/±45) m ] s lami- 
nates than for the [(±45/0/90) m ] s laminates. 

Unbalanced. Symmetrically Laminated Cylinders 

Unbalanced, symmetric laminates exhibit anisotro- 
py in the form of extensional-shear coupling ( A l6 * A 26 
* 0) in addition to flexural anisotropy (D, 6 * D 26 * 0 ). 
For these laminates, the value of the anisotropy parame- 
ter A given by Eqs. (20b) and (33) is not equal to unity. 


Simplifcd expressions for the anisotropy parameter sf {) 
and the first-order correction factor C v defined by Eqs. 
(33) - (35), arc given by Eqs. (C23) and (C24), respec- 
tively. Equation (C24) indicates that the value de- 
pends on coupling between the membrane and flexural 
anisotropies. 

Values of the orthotropy parameter 0 for [(+<|>) :m ] s 
symmetric, unidirectional laminates for the nine material 
systems considered herein are also presented in Fig. 4; 
that is, the curves presented in Fig. 4 for the [(±(j>) m ] s sym- 
metric angle-ply laminates are identical to those for the 
corresponding [(+<|>) 2 m ] s symmetric, unidirectional lami- 
nates. Thus, the orthotropy behavioral trends for the un- 
directional laminates are identical to those discussed 
previously for the symmetric angle-ply laminates, and 
are also independent of the stacking sequence number 
m. 

Results for the anisotropy parameter and the 
first-order correction factor are shown in Figs. 8 and 
9, respectively, for the [(+<J>) 2 J S symmetric, unidirection- 
al laminates with the nine material systems considered 
herein and are independent of the stacking sequence 
number m. The results in Fig. 8 show a substantial vari- 
ation in with fiber orientation and with material sys- 
tem. The results show that /^ 0 is the most pronounced 
for values of the fiber angle § between approximately 
55 deg and 80 deg, and that the contribution of the anisot- 
ropy to the attenuation behavior is essentially insignifi- 
cant (less than 1.05) for values of $ < 25 deg and <f> > 85 
deg. Moreover, the largest variation in /^ 0 with fiber an- 
gle is exhibited by the laminates made of the P-100/3502 
pitch-epoxy material and the smallest variation is exhib- 
ited by the laminates made of boron-aluminum material. 
Values of grange from approximately 1.42 for the max- 
imum point on the curve for the P-100/3502 pitch-epoxy 
material to a value of 1 . 

The results shown in Fig. 9 for the first-order cor- 
rection factor for the [(+<S>) :m ] s symmetric, unidirec- 
tional laminates indicate a substantial relative variation 
in with fiber orientation and with material system, but 
all of the magnitudes of are less than approximately 
0.45. Moreover, the magnitude of C\ is less than approx- 
imately 0.2 for all of the materials except the P-100/3502 
pitch-epoxy material. For the upper bound of thinness of 
thin-shell theory, given by h/R = 1/20, the contribution of 
to the anisotropy factor defined by Eq. (33) is practi- 
cally negligible. Equation (C24) indicates that the insig- 
nificance of S\ means that the coupling of the membrane 
and flexural anisotropies are negligible for these lami- 
nates. The insignificance of <?, is illustrated and verified 
in Fig. 10 for the [(H-<J>) 2m ] s symmetric, unidirectional 
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laminates made of IM7/5260 graphitc-bismaleimide ma- 
terial (black curves) and made of P-100/3502 pitch-ep- 
oxy material (gray curves), for h/R = 1/20. The finely 
dashed curves shown in Fig. 10 correspond to 90%-de- 
cay lengths for which the anisotropy is neglected. In 
contrast, the solid curves and the coarsely dashed gray 
curve include the effect of the membrane anisotropy and 
are shown for values of \l = 0, l, and 1.5. For these val- 
ues, results that correspond to the Sanders-Koiter theory 

3 

and the Love-Kirchhoff theory are given by ft = ^ and 

(a=l, respectively. Results that correspond to Donnell’s 
equations are given by ]X = 0. The solid curves in Fig. 
10 for p. = 1 and 1.5 are based on the exact solution that 
uses Eq. (20b) for the anisotropy factor. The correspond- 
ing curves that are based on the approximate formula for 
the anisotropy parameter that is given by Eq. (33) are 
identical. The solid curves and the coarsely dashed gray 
curve indicate that varying |i yields a small effect, which 
implies that all three shell theories yield essentially the 
same results and that /? 0 for the [(+<t>) 2 m ] s symmet- 
ric, unidirectional laminates. Comparing the solid and 
finely dashed curves in Fig. 10 also indicates that ne- 
glecting the membrane anisotropy underestimates the 
bending boundary-layer decay length, by as much as ap- 
proximately 3 1 % and 2 1 % for shell walls made of P- 100/ 
3502 pitch-epoxy and IM7/5260 graphite-bismaleimide 
materials, respectively. 

Values of the orthotropy parameter 0 for [(+45 2 /0/ 
90)J S and [(0/90/+45 2 )J s laminates made of the nine 
material systems considered herein are also presented in 
Fig. 7. More specifically, the values of 0 for these lam- 
inates are identical to the values for the corresponding 
quasi-isotropic laminates. Results for the anisotropy pa- 
rameter^ defined by Eq. (34) are shown in Fig. 1 1 for 
[(+452/0/90)^ and [(0/90/+45 2 )J s laminates made of 
the nine material systems considered herein. The results 
in Fig. 11 show no significant variation in /4 0 with the 
stacking sequence number m, and only a slight variation 
(less than approximately 9%) with material system. Val- 
ues of range between approximately 1.1 and 1. Cor- 
responding results for the first-order correction factor 
defined by Eq. (35), that are not shown herein, were ob- 
tained that indicate that all of the values of C\ for the 
[(+452/0/90) m ] s and [(0/90/+45 2 )J s laminates are less 
than approximately 0.1. These values indicate that the 
contribution of to the anisotropy factor defined by Eq. 
(33) is practically negligible. Thus, /t * for these 
laminates. The values of shown in Fig. 1 1 suggest that 
neglecting the anisotropy would, at most, underestimate 
the bending boundary-layer decay length by approxi- 
mately a 10%. The insignificance of C\ also means that 


the coupling of the membrane and flexural anisotropies 
are unimportant with regards to the primary effect of the 
individual shell anisotropies that is captured by the pa- 
rameter /^ n . 

Balanced. Unsymmetricallv Laminated Cylinders 

Balanced, unsymmetric laminates may, in general, 
exhibit anisotropy in the form of coupling between pure 
bending and twisting (D 16 * D, 6 * 0 ) and coupling be- 
tween membrane and bending action, which is manifest- 
ed by nonzero values for any of the subscripted B-terms 
in Eq. (A 15). These laminates do not, however, exhibit 
extensional-shear coupling ( A 16 = A 26 = 0). For the un- 
symmetric laminates that are discussed subsequently, the 
first ply in the stacking sequence is the innermost ply of 
a cylinder. Simplifed expressions for the anisotropy pa- 
rameter and the first-order correction factor <f t , de- 
fined by Eqs. (33) - (35), are given by Eqs. (C27) and 
(C29), respectively. Equations (C28) and (C29) indicate 
that the value depends on coupling between the flex- 
ural anisotropy and the anisotropy caused by unsymmet- 
ric lamination. 

Results for regular, antisymmetric angle-ply lami- 
nates are shown in Figs. 4 and 12-16. In particular, val- 
ues of the orthotropy parameter 0 for [(±<))) m ] T 
unsymmetric laminates made of the nine material sys- 
tems considered herein are also presented in Fig. 4; that 
is, the orthotropy-parameter curves presented in Fig. 4 
for the [(±<|)) m ] s symmetric angle-ply laminates are also 
identical to those for [(±4>)J T unsymmetric laminates. 
Thus, the orthotropy behavioral trends for the [(±<f>) m ] T 
unsymmetric laminates are identical to those discussed 
previously for the corresponding symmetric angle-ply 
laminates, and are also independent of the stacking 
sequence number m. 

Results for the anisotropy parameter/^ defined by 
Eq. (C27) are shown in Fig. 12 for 2-ply [±<(>] x unsym- 
metric laminates made of the nine material systems con- 
sidered herein. The results in Fig. 12 show a substantial 
variation in /4 0 with fiber orientation and with material 
system, and show that/^ 0 is the most pronounced for val- 
ues of the fiber angle $ between approximately 15 deg 
and 60 deg. Moreover, the largest variation in s? 0 with fi- 
ber angle is exhibited by the laminates made of the P- 
100/3502 pitch-epoxy material and the smallest variation 
is exhibited by the laminates made of boron-aluminum 
material. Values of /? {) range from approximately 0.75 
for the minimum point on the curve for the P-100/3502 
pitch-epoxy material to a value of 1. The results in Fig. 
13 show the variation in with the fiber angle 4> and 
the stacking sequence number m for [(±<J>) m ] T unsym- 
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metric laminates made of the P- 1 00/3502 pitch-epoxy 
material. These results show a rapid decline in the im- 
portance of sf {) > that is manifested by the curve moving 
closer to /f {) = 1, as the stacking sequence number in- 
creases. For m = 2, 0.95 1 . 

Results for the first-order correction factor are 
shown in Fig. 14 for 2-ply [±<|>] T unsymmetric laminates 
made of the nine material systems considered herein. 
The results in Fig. 14 also show a substantial variation in 
with fiber orientation and with material system. How- 
ever, the maximum magnitude of C\ is less than 0.07 for 
all of the material systems. Results are presented in Fig, 
15 that show the variation in g x with the fiber angle <j) 
and the stacking sequence number m for [(±4>) m ] T un- 
symmetric laminates made of the P-100/3502 pitch-ep- 
oxy material. These results show significant reductions 
in with an increase in the stacking sequence number. 

Overall, the results in Figs. 14 and 15 indicate that 
the contribution of to the anisotropy factor defined by 
Eq. (33) is negligible for the upper bound of thinness giv- 
en by h/R = 1/20, which means that /l ~ s? Q . Thus, the 
results in Fig. 12 for the two-ply [±<J>] T unsymmetric 
laminates indicate that neglecting the shell anisotropy 
overestimates the bending boundary layer, by as much as 
approximately 33% and 22% for shell walls made of P- 
100/3502 pitch-epoxy and IM7/5260 graphite-bismale- 
imide materials, respectively. The insignificance of 
also means that the coupling of the flexural anisotropy 
and the anisotropy caused by unsymmetric lamination is 
unimportant with regards to the primary effect of the in- 
dividual shell anisotropies. The insignificance of is il- 
lustrated in Fig. 16 by the gray and by the black curves 
for the laminates made of P-100/3502 pitch-epoxy and 
IM7/5260 graphite-bismaleimide materials, respective- 
ly. The solid black and gray curves are for the upper 
bound of thin-shell theory that is given by h/R = 1/20. 
The finely dashed curves shown in the Fig. 16 corre- 
spond to 90%-decay lengths for which the anisotropy is 
neglected. In contrast, the solid curves include the effect 
of the shell anisotropy and are shown for values of |i = 0, 
1, and 1.5. The solid curves for ji = 1 and 1.5 are based 
on the exact solution that uses Eq. (20b). The corre- 
sponding curves that are based on the approximate for- 
mula for the anisotropy parameter that is given by Eq. 
(33) are identical. The solid curves indicate no signifi- 
cant effect of varying fl, which implies that all three shell 
theories yield essentially the same results for the [±<f>] T 
unsymmetric laminates. For [(±<()) m ] T unsymmetric lam- 
inates with m > 1 and made from any of the nine material 
systems considered herein, the results in Figs. 12 through 
15 indicate that neglecting the shell- wall anisotropy will 
have a small effect on the calculation of the bending 


boundary-layer decay length. 

Values of the orthotropy parameter 6 and the 
anisotropy parameter for (0 p /90 t| ) r unsymmetric 
cross-ply laminates are shown in Figs. 17 and 18 for the 
nine material systems considered herein and as a func- 
tion of the percentage of zero-degree plies. For this class 
of laminates, Eq. (20b) simplifies to Eq. (34); that is, a f- 
/4 ky This simplification means that the anisotropy param- 
eter is independent of p, which means that all three of the 
shell theories considered herein yield identical results. 

The results in Fig. 17 show a large variation in O 
with the percentage of zero-deg plies for most of the ma- 
terial systems. In addition, the results show a large vari- 
ation in 0 with material system for the laminates that are 
dominated by ninety-deg plies (less than approximately 
10% zero-deg plies) and by zero-deg plies ( more than 
approximately 80% zero-deg plies). Values of 0 vary the 
most for laminates made of P-100/3502 pitch-epoxy ma- 
terial, with values that range from approximately 0.3 to 
2.93. Most of the materials exhibit values of O in the 
range of approximately 0.5 to 2. 1 . 

The results in Fig. 18 also show a large variation in 
with the percentage of zero-deg plies for most of the 
material systems, and show a large variation with mate- 
rial system for laminates with less than 70% zero-deg 
plies. Moreover, the results show that/^ is the most pro- 
nounced (most different from a value of 1) for laminates 
with approximately 15% to 30% zero-deg plies. The 
largest variation in/^ 0 with percentage of zero-deg plies 
is exhibited by the laminates made of the P-100/3502 
pitch-epoxy material and the smallest variation is exhib- 
ited by the laminates made of boron-aluminum material. 
Values of /grange from approximately 0.57 for the min- 
imum point on the curve for the P-100/3502 pitch-epoxy 
material to a value of 1. Thus, in some cases neglecting 
the shell wall anisotropy overestimates the bending 
boundary layer, by as much as approximately 75% for a 
shell wall made of P-100/3502 pitch-epoxy material. 
This result is illustrated in Fig. 19 by the gray curves. 
Similar results are presented in Fig. 19 for (0 p /90 q ) T un- 
symmetric cross-ply laminates made of IM7/5260 graph- 
ite-bismaleimide material (black curves). The solid 
black and gray curves include the effect of the shell 
anisotropy and the Finely dashed curves shown in the fig- 
ure correspond to 90%-decay lengths for which the 
anisotropy is neglected. The results in Fig. 19 show that 
including the effect of anisotropy is particularly impor- 
tant for laminates with less than approximately 70% 
zero-deg plies. 

Unbalanced. Unsymmetricallv Laminated Cylinders 

Unbalanced, unsymmetric laminates may, in gen- 
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cral, exhibit full anisotropy in the form of coupling be- 
tween pure bending and twisting (D lfi * D 26 * 0 ) and 
coupling between membrane and bending action, which 
is manifested by nonzero values for any of the subscript- 
ed B-terms in Eq. (A 15), and extensional-shear coupling 
( A lfi * A 26 * 0). The expressions for the anisotropy pa- 
rameter/^, and the first-order correction factor & x that are 
given by Eqs. (C2) - (C22) indicate that/^, exhibits cou- 
pling between the membrane anisotropy and the anisot- 
ropy that is caused by unsymmetric lamination, and that 
& x exhibits coupling between all three types of anisotro- 
pies. One family of laminates that exhibits all of these 
anisotropies is the (70 p /0 q ) T unbalanced, unsymmetric 
laminates with p > 0 and q * 0. 

Values of the orthotropy parameter 0 and the 
anisotropy parameter s4 0 for (70 p /0 q ) T unbalanced, un- 
symmetric laminates are shown in Figs. 20 and 21, re- 
spectively, for the nine material systems considered 
herein and as a function of the percentage of seventy-deg 
plies. The results in Fig. 20 show a large variation in 0 
with the percentage of seventy-deg plies for most of the 
material systems. The results also show a large variation 
in 0 with material system for the laminates that are dom- 
inated by zero-deg plies (less than approximately 20% 
seventy-deg plies). Values of 0 vary the most for the 
laminates made of P-100/3502 pitch-epoxy material, 
with values that range from approximately 0.5 to 3. 

The results in Fig. 21 also show a substantial varia- 
tion in /4 0 with the percentage of seventy-deg plies for 
most of the material systems, and a large variation with 
material system for laminates with between approxi- 
mately 45% and 100% seventy-deg plies. The largest 
overall variation in /f 0 with percentage of seventy-deg 
plies is exhibited by the laminates made of the P-100/ 
3502 pitch-epoxy material and the smallest variation is 
exhibited by the laminates made of boron-aluminum ma- 
terial. Values of /4 0 range from approximately 1.4 to 
0.95, which correspond to the maximum and minimum 
points, respectively, on the curve for the P-100/3502 
pitch-epoxy material. 

Results for the First-order correction factor were 
also obtained for (70 p /0 ) T unbalanced, unsymmetric 
laminates made of the nine material systems considered 
herein, but are not included in the present paper. These 
results also show a substantial, relative variation in @ x 
with the percentage of seventy-deg plies, but overall the 
magnitude of C x is less than approximately 0.25 for the 
P-100/3502 pitch-epoxy material and less than 0.1 for 
the other materials. These results indicate that the con- 
tribution of & x to the anisotropy factor defined by Eq. 
(33) is negligible for the upper bound of thin-shell theory 


that is given by h/R = 1/20, which means that rt - zt {Y 
Thus, the results in Fig. 2 1 suggest that in some cases ne- 
glecting the shell-wall anisotropy may overestimate the 
bending boundary-layer decay length and in other cases, 
may underestimate the decay length. The insignificance 
of d x also means that the contribution of the flexural 
anisotropy to the coupling of the anisotropies is negligi- 
ble. The insignificance of & x is clarified in Fig. 22 for 
laminates made of P-100/3502 pitch-epoxy material 
(gray curves) and of IM7/5260 graphite-bismaleimide 
material (black curves). The solid black and gray curves 
are for the upper bound of thinness given by h/R = 1/20. 
The finely dashed curves shown in the figure correspond 
to 90%-decay lengths for which the anisotropy is ne- 
glected. In contrast, the solid curves include the effect of 
the shell anisotropy and are shown for values of [X = 0, 1 , 
and 1 .5. Moreover, the solid curves for p = I and 1 .5 are 
based on the exact solution that uses Eq. (20b). The cor- 
responding curves that are based on the approximate for- 
mula for the anisotropy parameter that is given by Eq. 
(33) are identical. The solid curves indicate a negligible 
effect of varying p, which verifies that and im- 

plies that all three shell theories yield essentially the 
same results for the (70 p /0 q ) T unbalanced, unsymmetric 
laminates. In addition, the results show that neglecting 
the shell wall anisotropy, for the most part, underesti- 
mates the bending boundary -layer decay length, by as 
much as approximately 1 6% and 6% for shell walls made 
of P-100/3502 pitch-epoxy and IM7/5260 graphite-bis- 
maleimide materials, respectively, and with approxi- 
mately 20% seventy-deg plies. In addition, the results in 
Fig. 22 show that neglecting the shell wall anisotropy un- 
derestimates the bending boundary-layer decay length, 
by as much as approximately 31% and 20% for shell 
walls made of P-100/3502 pitch-epoxy and IM7/5260 
graphite-bismaleimide materials, respectively, and with 
approximately 100% seventy-deg plies. There is only a 
very small range shown in Fig. 22 where neglecting the 
shell wall anisotropy overestimates the bending bound- 
ary-layer decay length, and for this region, the effect is 
negligible. 

Concluding Remarks 

An analytical study of the attenuation of bending 
boundary layers in both balanced and unbalanced, sym- 
metrically and unsymmetrically laminated-composite, 
thin cylindrical shells has been presented for nine con- 
temporary material systems. The analysis is based on the 
linear Sanders-Koiter shell equations and contains the 
Love-Kirchhoff shell equations and Donnell’s equations 
as special cases. With this analysis, two nondimensional 
parameters have been indentified that characterize and 
quantify the effects of laminate orthotropy and laminate 
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anisotropy on the bending boundary-layer decay length 
in a very general and encompassing manner. The anisot- 
ropy parameter includes the effects of anisotropy in the 
form of coupling between pure bending and twisting that 
appears in many symmetric laminates to some extent, 
coupling between extension and shear that is present in 
unbalanced laminates, and coupling between membrane 
and bending action that is present in unsymmetric lami- 
nates, 

A substantial number of structural design technolo- 
gy results for the bending boundary-layer decay length 
have been presented for a wide range of laminated-com- 
posite shell structures that should be useful additions to 
the structural designer’s collection of preliminary design 
tools. Moreover, the analysis and results should provide 
additional physical insight into the fundamental behavior 
of general laminated composite shell structures and pro- 
vide a common basis for assessing bending boundary- 
layer attenuation for the vast range of laminate construc- 
tions that are possible. Furthermore, the results should 
be useful for the design of specimens for material char- 
acterization tests, for instrumenting structural verifica- 
tion tests, and for defining finite-element meshes. For ail 
the laminate constructions considered in the present 
study, the results show that the differences between re- 
sults that were obtained with the Sanders-Koiter shell 
equations, the Love-Kirchhoff shell equations, and Don- 
nell’s equations are negligible. The results also show 
that the effect of anisotropy in the form of coupling be- 
tween pure bending and twisting has a neglible effect on 
the size of the bending boundary-layer attenuation length 
of the balanced, symmetrically laminated cylinders con- 
sidered. Moreover, the results show that the coupling of 
the membrane and flexural anisotropy and the anisotropy 
caused by unsymmetric lamination is generally unimpor- 
tant with regards to the primary effect of the individual 
shell anisotropies on the bending boundary-layer decay 
length. The only exception encountered was for unbal- 
anced, unsymmetrically laminated cylinders for which 
coupling of the membrane anisotropy and the anisotropy 
caused by unsymmetric lamination is a primary effect, as 
expected. The results also show that in some cases ne- 
glecting the shell anisotropy results in underestimating 
the bending boundary-layer decay length and in other 
cases it results in an overestimation. 

Appendix A: Sanders-Koiter Eq u ations 

The linear Sanders-Koiter shell equations* 6,17 are 
presented in this appendix for a right-circular cylinder 
with a radius that is given by R. For these equations, x 
and 9 denote the axial and circumferential coordinates, 
respectively. First, the equilibrium equations are pre- 
sented, then the kinematic equations and the constitutive 


equations arc presented. Last, the boundary conditions 
are given for a complete right-circular cylinder at an 
edge that is given by a constant value of the axial coordi- 
nate, x. 

Equilibrium Equations 

The equilibrium equations are given in a form sim- 
ilar to those found in Ref. 22; that is, 


dN, 1 3N,„ c, 3M, 9 n 

~3T + R~3e"~2R J “59“ q ' 

(Al) 

1 3N„ c,_ . c 2 3M, 9 _n 

~5T + R“3r + R Qb + 2R1T + q * “ 0 

(A 2) 

3Q, 1 3Q„ N„ _ q 

"35T + IC3e R +q "- u 

(A3) 

3M, j 3M„, /-) _ q 

^r + R“5e" Q ‘"° 

(A4) 

3M i 3M 9 p, _ q 

ut + r\w Q »-° 

(A5) 


where N x , N e , and N x0 are the membrane stress result- 
ants; Q x and Q 0 are the transverse shear-stress resultants; 
M x , M e , and M x0 are the bending stress resultants; q x , q 0 , 
and q n are the applied surface tractions; and c, and c 2 
are constants that identify the equations of other shell 
theories that are considered herein. In particular, the 
Sanders-Koiter equations are given by c } = c 2 = 1 and 
the Love-Kirchhoff equations are give by c, = 1 and c 2 = 
0. Donnell’s equations are given by c t = c 2 = 0. This 
convention is used throughout the present study. 

Kinematic Equations 

The kinematic equations are given by 


e o 9u 

(A6) 

pH _ 1 9v ,w 
fca “ R30 T* 

(A7) 

_ 3v 1 3u 
~ 3x R3e 

(A8) 

R» _ 3w 

P-- _ W 

(A9) 

B " = c, 1 3w 

Pe R R7)6 

(A 10) 

r° _ c :/ 3v 1 3u ) 

p "-Tl3I"R30j 

(All) 
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3P" a 
= u^~~ 


3x J 


v-.. _ ± 3 P« _ £i3v _ I 3 2 w 

k “ _ R^ 0 R : 3e R : 30 2 


(A 1 2) 
(A 13) 


where u(0), v(0), and w(0) are applied edge displace- 
ments; P(9) is an applied edge rotation; and N s (9), 
T(9),V(0), and M„(9) are applied edge loads. 




ap; 
+ ~zc 


2 a 2 w 

R5x59 


4' 


(c, + )c,’ 


3v 

5x 


c. 


3u 

2R : ae 


(A 14) 


where u, v, and w are the axial, circumferential, and nor- 
mal displacements of a point of the shell middle surface; 

£°, and 7!» are the membrane strains; P°*, Pa, and 

Pn are the rotations; and Kj t and Ke are the bend- 
ing strains. The displacement w is positive when it is 
outward from the cylinder reference surface. 


The isothermal constitutive equations are given in 
matrix form by 


Appendix B: Equations for Axisymmetrv 

The linear Sanders-Koiter shell equations that are 
presented in Appendix A for a right-circular cylinder 
with a radius R are specialized in this appendix for the 
case of axisymmetric behavior. For these equations, x 
and 0 denote the axial and circumferential coordinates, 
respectively. The specialization to axial symmetry is 
conducted by eliminating all terms in the equations of 
Appendix A that are differentiatied with respect to the 
circumferential coordinate, 0. First, the equilibrium 
equations, the kinematic equations, and the constitutive 
equations are presented. Then, the boundary conditions 
are given for a complete right-circular cylinder at an 
edge that is given by a constant value of the axial coordi- 
nate, x. Last, the axisymmetric equations are manipulat- 
ed into a single ordinary differential equation in terms of 
the normal displacement w(x). 



A„ 

a 12 

A„ 

B„ 

B« 

b, 6 

a b 

A« 

a 26 

B,: 

b 22 

b 26 

a 16 

A* 

A„ 

b 16 

B;„ 

B* 

B i, 

B« 

b, 6 

pn 

d, 2 

d, 6 

B I2 

b 22 

b 26 

Bb 

d 22 

d :6 

B* 

b 2# 

B« 

d, 6 

d 26 

D* 



(A 15) 


where the subscripted A, B, and D terms of the matrix 
are the stiffnesses of laminated composite shells that are 
obtained from the Love-Kirchhoff shell theory. More- 
over, the constitutive terms in Eq. (A 15) are identical to 
those for laminated-composite plates that are given in 
Ref. 18, p. 198. 


Equilibrium Equations 

The equilibrium equations for axisymmetric behav- 
ior are given by 


dN x , . n 

"dx 1 + q* x = 0 


dN„ e c 2 dM ^ + ^ \ q 

dx R^ 9 2R dx 


q n(x)=0 


dM 

dx 


-~Q* = o 


(Bl) 

(B2) 

(B3) 

(B4) 


Boundary Conditions 

The boundary conditions for an edge that is defined 
by a constant value of the axial coordinate x are given by 


N, = N,(0) or 

u = 0(0) 

(A 16) 

N ’ ,b + r( C| + = ^(9) 

or v = v(0) 

(A 17) 


or w = w(0) 

(A 18) 

M, = R(0) or 

P°. = P(0) 

(A 19) 


^jf-Q0 = O (B5) 

where the membrane stress resultants N x , N 0 , and N x9 ; 
the transverse shear-stress resultants Q x and Q e ; the 
bending stress resultants M x , M e , and M xe ; and the 
applied surface tractions q x , q e , and q n are functions of 
only the axial coordinate, x. 



The kinematic equations are given by 
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(B 17) 


r< > _ du 
* dx 

(B6) 

e o = w 
fc ° R 

(B7) 

y 1 -dx 
y * dx 

(B8) 


(B9) 


(BIO) 

no _ C 2 (Jv 

p " _ 2dx 

(Bll) 

K° = — = - d ' W 

* dx dx 2 

(B12) 

K°e = 0 

(B 13) 

K „ _i S o + i c 

k .»- r P„+ dx - rY |+ 2° 2 

)e <® 14 > 

where the middle-surface displacements u, v, and w; the 
membrane strains £° , £ 0 , and yl * ; the rotations P° x , P° 0 , 

and P° ; and the bending strains k° and k ° 6 are func- 
tions of only the axial coordinate, x. 


N * + r( C, + 2 C; ) M '“ = T or v = v 

Q x = V or w = (B 18) 

M, = M, or (3>P (B 19) 


where the applied edge displacements 0, v, and w; the 
applied edge rotation P ; and the applied edge loads N s , 
T, V, and M, are all constants. 


The bending boundary-layer equation is obtained 
by first noting that integration of Eq. (Bl) yields 


N s = - J q x dx + C = N(x) 


(B20) 


where C is a constant of integration that is determined 
from the boundary condition given by Eq. (B16). Next, 
Eqs. (B2) and (B5) are combined to get 


dN xa 

dx 





dM 

dx 


x6 


+ q fl (x) = 0 


(B21) 


For convenience, the parameter 


p. = Ci + ^c 2 (B22) 


Constitutive Equations 

The isothermal constitutive equations reduce to 



A,, 

a i2 

a, 6 

B„ 

b, 2 

B 

A„ 

a 22 

a* 

b 12 

b 22 

B~ 

a 16 

A,, 

a„ 

b 16 

b* 

B, 

B„ 

b, 2 

b, 6 

D„ 

d, 2 

D 

b 12 

b 22 

B* 

d 12 

D n 

D 

b 16 

b 26 

B« 

D ls 

D* 

D ( 



(B15) 


is introduced such that the Sanders-Koiter equations are 

'i 

given by = ^ and the Love-Kirchhoff equations are 

given by \1 = 1. Donnell's equations are given by \x = 0. 
Similarly, the function 

T(x) = N i9 + £m (6 (B23) 

is introduced so that Eq. (B21) becomes 


where the subscripted A, B, and D terms of the matrix 
are the usual constitutive terms of classical Love-Kirch- 
hoff-type laminated composite shell theory or classical 
laminated plate theory (e.g., see p. 198 of Ref. 18). 


g + q 6 (x) = 0 (B24) 

and the corresponding boundary condition given by Eq. 
(B17) becomes 


The boundary conditions for an edge that is defined 
by a constant value of the axial coordinate x are given by 


T = T or v = v 
Integration of Eq. (B24) yields 


N, = N S or u = u (B16) 


T = - J q a dx + C = T(x) 


(B25) 


(B26) 


15 

American Institute of Aeronautics and Astronautics 


where C is a constant of integration that is determined 
from the boundary condition given by Eq. (B25). Next, 
Eqs. (B3) and (B4) are combined to give 

^-^ + q„(x) = 0 (B27) 

The next step in the analysis is the simplification of 
the constitutive equations. First, by using Eqs. (B8) and 
(B22), Eq. (B 14) is expressed as 


and Eq. (B26) is expressed as 

S|,l dx +S “¥ + S “dx _B, ',iF“T (x) = 0 (B37) 


du 


Equations (B36) and (B37) are then solved for ^ and 


dv 

dx 


to get 


C = (B28) 

By using Eqs. (B23) and (B28), the constitutive equa- 
tions are expressed as 


du _ S M M(x)-S, t T(x) 
dx a„S„-S' 6 

+ —= — T> t 13 - 38 ) 


and 


N . 

, 

An A 12 A, 6 B 

N e ' 

1 _ 

A 12 A 22 A 26 B 

T 1 

f ~ 

a 16 a 26 a^b 

M.J 

1 

®11 ®12 ®16 ^ 



(B29) 


dv _ A„T(x)-A, 1 N(x) 
dx a 


(A l; S, 6 -A„S a )f + (a ii B 14 -A,3h) 


d~w 

dx J 


11^66 16 


(B39) 


® 12^x 4" ^22^0 + ^26^6 + Dl2^x 

(B30) 

where 


Xi6 = A i« + H(g;)-ff 

(B31) 

x -a +ufh.]MM 
26 — 26 ^ El 1 

(B32) 

S» = A < . + 2n(£)5 f . ll ’(i) I £f 

(B33) 

b >= b “ + Ke)t 

(B34) 


(B35) 

The motivation for writing the constitutive equations in 
this form is that the matrix equation given by Eq. (B29) 
is the only part of the of the full constitutive equations 
that appear in the strain energy density function, which 


is used in the present paper to determine the correspond- 
ing positive-definiteness conditions. With these simpli- 
fied constitutive equations and Eqs. (B6) - (B8) and 
(B12), Eq. (B20) is expressed as 


Equation (B39) indicates that the circumferential dis- 
placement v(x) becomes uncoupled from the axial dis- 
placement u(x) and the normal displacement w(x) 

when S t6 = A 26 = B l6 = 0, which implies that A 16 = A 26 
- B 16 = B 26 = D 16 = 0. In addition, the constitutive equa- 
tion, Eq. (B29), indicates that N x , N e , and M x become 

uncoupled from the torsional, shear strain Yxe when 
A i6 = A 26 = B 16 = 0 , and that T, that is defined by Eq. 

(B23), becomes uncoupled from £ e> and K. Fur- 
thermore, Eq. (B30) indicates that M 0 becomes uncou- 
pled from Yq when B 26 = 0 , which implies B 26 = D 26 = 
0. 

Next, Eqs. (B38) and (B39) are then substituted 
into Eqs. (B6) and (B8), and the resulting expressions for 

and along with Eqs. (B7) and (B 1 2) are substitut- 
ed into the constitutive equations, (B29). This action 
converts the strains and stress resultants in Eq. (B29) into 
functions of the displacement w(x). Substituting the ex- 
pressions for N e and M x into Eq. (B27) yields the bend- 
ing boundary-layer equation that is given by 

c, $ +c 4f +C)W=C 4(x) (B40) 


Ai >d^ +Ai 3 f +A -S' B "I ? _R(X)=0 (B36) 


The constant coefficients are given by 
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— D n 


1 - 


(a u A„-A’,)d„ 


(B41) 


C, = - wB 


2- R "I2 

/ 


|A 1a A » — A i 2 A w |b j] + | A )2 A , A A,, A 


A„A„- A ! I5 


(B42) 


(a,,Ai, - A ujAjs A„A m A Ii A 1 .+ 2A, ; A^ JS 
C 3 = ' t pr (B43) 

r j (a u a„- a’ 6 ] 

The function C 4 (x) is given by 
C 4 (x) = q„(x) + 

( a lt A ,6 - A i2 A m ]N(x) + ( A „A - A „ A 26 ]T(x) 

r(a„a«-a; 6 ) 

[b,iA m - A 16 B l6 ) | A 2 |B 16 — B,,A i A j 


(a„A„- A‘, 6 ) 


(B44) 


These expressions are simplifed further by introducing 
the following expressions 


a„ = 


a 22 = 


a* = 


A ia A m A 12 A w 


(A„A a A^Ja^ A,,A 26 A 22 A I6 + 2 A l; A i4 A ^ 

AnA^-A', 

[A„A22- A“ 2 |A66- A n A^6 — A 2 2A’ 4 + 2A,;A,6A26 

A | 2 A J6 — A ,,A 26 


(B45) 


(B46) 


(B47) 


(a „A 2 j - A^A* - A ..A* - AjjA 1 ,. + 2A , ; A 16 A 

6 2 , = — |a i; B n + SjtBu + SisBuj (B48) 

, A^;, + A ll Bj 6 -2A le B, 1 B l6 ,„ jr>N 

e = 1 7 — Z7\ (B49) 

(a„a 66 -S I4 Jd 11 

By using Eqs. (B45) - (B49), Eqs. (B41) - (B43) are 
expressed as 

C, =D n e 


P — 2_ ^:i 

J 'Ra„ 


(B50) 

(B51) 


r = ! 

Vi n 2, 


R~a Yl 


(B52) 


Similarly, For the case where the second derivatives of 
N(x) and T(x) are zero valued, Eq. (B44) becomes 


C 4 (x) = q n (x) + 


a, ; R(x) + a,J(x) 

Ra 22 


(B53) 


The desired form of the bending boundary-layer equa- 
tion is obtained by dividing Eq. (B40) by C t ; that is, 


£ + 4 S 3 # t 4 < 5 w = p < x > 


.d 2 w 


(B54) 


where the constants S and Q are given by 

S = - 5 21 

4C, 2Ra l2 D n e 


Q = — = 


i 


(B55) 


(B56) 


4C, 4R 2 a 22 D,,e 

The function P(x) is given by 

P (JC) = c «w = qn(x) , a, 2 N(x) + a 26 T(x) (B57) 

DuC ^a 22 D j 

for the special case when the second derivatives of 

N(x) and T(x) are zero valued. The quantity D n e that 
appears in Eqs. (B54) - (B56) is sometimes referred to, 
in some contexts, as a reduced bending stiffness. 12 

A ppendix C: Anisotropy-Factor Eq uations 

The first-order approximation of the anisotropy 
factor/# that is used herein is given by 




\+\ie 


trJj 


(Cl) 


where /# 0 is the value of Eq. (20b) with p = 0 in the 
terms with the overbars. This expression is given by 


^0 = 


(A u A 2? -A 2 2 ) 


^22^n 


1 - 


b 2 j 


\J a^Dn c Q 


(C2) 


which is the anisotropy factor that corresponds to Don- 
nell’s equations. The terms a 22 , b 2l , and e 0 are given by 
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A l6 A i6 - A.jA^ 

(a„A„ - A',) 

Aiy>™ A 1 1 A 2ft — A 32 A , A +■ 2 A ti A l6 A 2A 
A A - A 1 

^ V, •* !6 //^ A\ 

(a.^-a 1 ,,) 

) 

A.S6 — A U A 2A — A 22 A , 6 + 2 A tJ A lft A 2A 
A ti A ia — A j j A 26 

(A„A n - A’j) 

b„=- 

(C5) 

A m - A„A :6 - A„A ,. + 2A 1 j A, 4 A j , 

( a i2®ll ^22® 12 ^2S® l<v) (C6) 


f’is “ A l6 | A :*A \ i A ia] “b 

A i^A. n A 22 + A 12 j — 2A„A 1: A 2t 

= 2^A i2 A 16 — A n A 26 j^A i^A 22 ~ A l2 A 26 j 
§J6 = (a |2 A 16 — A II A 26 A |<yA > 6 “ A ^A^j 
£?26 == (A 12 A | 6 — A nA^ftj^A I, A 66 — A , 6 j 


(Cl 5) 

(Cl 6) 
(Cl 7) 
(Cl 8) 


A^;, + A u B| 6 -2A 16 B u B 16 
(a ii A 86 — A , 6 jD n 


(C7) 


The term is a first-order correction to the results that 
correspond to Donnell’s equations and is given by 


V a 2iDii*o a 2:*i + <Wn) + 2a 22 ( b J2 ^ 0 b 21 ^| 2a| 21 bj t £ 0 

e ,= r 1 1 (C8) 

4a 22 ^ 0 \J a 22 D,,* 0 — b 2 , 


where 

- ^(^ 16 ^ n ~ 


C\ = 


D 16^^ 1 1 A 66 “ Ai 6 ) + 


(AhA^-A^J D u h 

“ B i6( B n A 66 + A M B 66 j + A 16 (a U B 66 + Bj 6 j 


3-i i? *“ 


mAjj- a! 2 ]a 


B lt/|6 + B 2^26 + B 


(C9) 

(CIO) 


Sft6 — [A u A| 6 A [ [ A 26^ (Cl 9) 

h 16 = A 12 A66^A l2 A66 — A | 2 j + A,|A26^A l2 A26- A^Av,) 

■b A 1& A 22 ^A l2 A 16 — A nA :6 j (C20) 

h 2 6 = A nA^A n A 2 2 — A j2 j + A ,6^ A | jAtj ~ 2 A , 2 j 

“b A j j A 26^2 A 12 A | 6 A | j A 26 | (C21) 

^66 = 2^A i j A 26 Aj 2 A i6 j^A nA 22 — A l2 j (C22) 


Special Cases for j4^ and ^ 

Simplifications to /4 0 and & K are presented below 
for unbalanced and balanced, symmetric laminates and 
for balanced, unsymmetric laminates, that include the 
subclasses of general antisymmetric laminates, antisym- 
metric cross-ply laminates, and antisymmetric angle-ply 
laminates. 


a m — 2 




(a itAjj- A’ 2 ]i 


I A 66 ~ A , | A jg — A 2jA + 2A l2 A )6 A 24 


(C 11) 




Bjfth i6 * B 2A fr 26 + 


( A i jA 22 — A^ 2 |a — A n Al 6 - A 23 Al 6 4-2A, 2 A| ft A 2 


(Cl 2) 


^ 121 — I a n 2 B u + a 122 B t2 + a, 2ft Bi 6 + ^ I (C13) 


Unbalanced and balanced, symmetric laminates. 
For unbalanced, symmetric laminates, A l6 ^ 0, A 26 * 0, 
and B n = B 12 = B 22 = B 16 = B 26 = B 66 = 0. For this special 
case, 


^0 = 


(A t ]A 2 2- Ai^A.jAj*- A , 6 j 
A jJ (a nAjj — A i^A^- A | |A 26 - A^A 16 + 2A t2 A |*A 2( 


(C23) 


and 

^t6 = A 26 ^A 16 A 22 — A n A 26 j + 


A 2 6^ A n A 22 + A , 2 j 2A 12 A 22 A,6 


(Cl 4) 


which agrees with the corresponding equations given by 
Reuter 4 , and 


e x 


ifef A [2 A |6 A u A 26 j a 

P (a„a 22 -ai 2 )/ 

Dh(AhA 66 - A] fi ) 
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- A hA^ - A 22^*16 + 2A 12 A lft A 


(C24) 


0 in addition to the shear-cxtensional coupling terms. 
For this special ease, where 


For balanced, symmetric laminates, A l6 = A 26 = 0 in 
addition to the subscripted B-matrix constitutive terms. 
For this special case, = 1 and S\ = 0. 

Balanced, unsvmmetric laminates. For balanced, 
unsymmetric laminates, A I6 = A 26 = 0, which yields the 
following simplified expressions 


_■ K 

0 A ,.D , 


(C33) 


which agrees with the corresponding equations given by 
Reuter 4 , and 


b;, 


b: 


e ° 1 A„D|, AJ3, 


A = 


A i iB |2 A 12 B n 


y A|,D n (A lt A 2 j A 12 )e 0 

^0 = + a] 


B, t 

AJi 


Dl6 ®11®I6 ® )6 B 66 


D„ A u D n A 66 D l 


“ 4(1+ A) 


— (1+2A) + 


2B, 6 ^A„B» A 12 B I6 J 

A u D n | A n A» ~ A uJ«o 


(C25) 

(C26) 

(C27) 

(C28) 

(C29) 


For the subclass of balanced, antisymmetric laminates, 
D 16 = D 26 = 0 in addition to the shear-extensional cou- 
pling terms, which yields the following simplification 


2B; 6 

B„ 

A*D„h 

A„ 


(C30) 



Further simplifications can be made to Eqs. (C31) 
and (C32) for [0/90/.. ./90] antisymmetric-cross-ply- 
laminate shell walls with an even number of layers that 
have identical material properties. For these laminates, 
the plies are specially orthotropic and their principal ma- 
terial directions are oriented at 0 deg and 90 deg to the 
cylinder axes in an alternating manner. In particular, the 
major principal axes of the odd-numbered and even- 
numbered plies are aligned with the x- and 0-axis, re- 
spectively, with the first ply in the stacking sequence lo- 
cated at the inner surface of the cylinder. Moreover, the 
odd-numbered plies have the same thickness and the 
even-numbered plies have the same thickness, but these 
two thicknesses are, in general, different. The laminate 
stiffnesses are given in Ref. 18 (see pp. 224-226) in terms 
of the number of layers N, the thickness ratio M, the ratio 

E 

of the principal elastic moduli F = g 1 , for which 0 < F 

< 1, and the reduced, plane-stress lamina siiffnesses. 
The thickness ratio is defined by 


that is applied to Eq. (C29). For the subclass of (bal- 
anced) antisymmetric cross-ply laminates, B l2 = B 16 = 
B 26 = B 66 = 0 , B 22 = -B u , and D I6 = D 26 = 0 in addition 
to the shear-extensional coupling terms. For this special 
case, = 0 and 


< C35 > 

where t (k) denotes the thickness of the k" 1 ply and 



A|;B,, 


\J AiiDh(A u A 2 2 A 12 )^ 0 


(C31) 


where 

, Bu 

'r'-IT (C32) 

For the subclass of balanced, antisymmetric angle- 
ply laminates, B n = B l2 = B 22 = B 66 = 0 and D l6 = D 26 = 


kk) ^(k) 2 (* m + 1<2) ) (C36) 

is the total laminate thickness. For the antisymmetric 
cross-ply laminates, t U) and t (2) are the thicknesses of the 
0-deg and 90-deg layers, respectively, and are denoted 
herein by ^ and t^, respectively. Substituting the non- 
zero laminate stiffness expressions for this class of anti- 
symmetric cross-ply laminates that are given in Ref. 18 
into Eqs. (21), (C31), and (C32) yields 
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c > = 


[ 1 - ( I - F)Q](M + 1 )(M + F) 


(M + F)(I +MF)-[(I + M)Fv i2 ]‘ 


(C37) 


• = V^oj 


1 + 


VTZMF(1-F)v, 2 


NI 


(C38a) 


where 


E= [1 -(1 -F)Q](M + I)(M + F)x 


{(M + F)(l + MF) - [(M + I)Fv 12 ] 2 }e 0 


(C38b) 


12M 2 (1-F) 2 

N 2 (M+ 1) 5 (M + F)[1 -(1 -F)Q] 


(C39) 


1 8M(M- 1) 

^ 1+M N 2 (M + l) 3 


(C40) 


and v, 2 is the major Poisson’s ratio. For the special, but 
practical, case of regular antisymmetric cross-ply lami- 
nation, all plies have the same thickness and Eqs. (C37) 
- (C39) reduce to 


0 = 


1 - 



(C41) 



1 + 


vTIv i2 F(1 -F) 


N(1 + F) v /[(l+F) : -(2Fv li ) ! ] e „ 


(C42) 



(C43) 
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Table 1: Lamina properties. 


r £:*• - - .. 

I Lamina 

. . . Material Systems ;; ; , 


S-glass- 
~ epoxy 

Kevlar 49- 
epoxy 

IM7/ 

5260 

AS4/ 

3502 

AS 4/ 
3501-6 

Boron- 

epoxy 

IM7/ 

PKTI-5 

P-100/ 
3502 . 

33 

7.5 

11.02 

22.1 

18.5 

20.01 

29.58 

20.35 

53.5 

21 

1.7 

0.8 

1.457 

1.64 

1.30 

2.68 

1.16 

0.73 

nSrn 023 

0.25 

0.34 

0.258 

0.30 

0.30 

0.23 

0.29 

0.31 


0.80 

0.33 

0.860 

0.87 

1.03 

0.81 

0.61 

0.76 

BSiEll 3 - 2 

3.5 

-2.22 

0.0125 

0.25 

-0.167 

3.38 

-0.14 

-0.64 

110 

11.0 

43.89 

14.91 

16.2 

15.6 

16.83 

16.85 

17.2 


* The subscripts I and 2 denote the longitudinal (fiber) and transverse (matrix) directions of a specially orthotropic lamina, respectively. 


Table 2: Comparison of results for specially orthotropic materials with Ref. 6. 



* The subscripts 1 and 2 denote the major and minor principal directions, respectively, of the specially orthotropic materials defined 
in Ref. 6. 

t The quantity IRcC/h} is defined by Eq. (37). 
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Fig. 1 Nondimensional 90% -decay length for symmetrically and unsymmetrically laminated cylinders, as a 
function of laminate orthotropy. 



Anisotropy parameter,/# 


Fig. 2 Nondimensional 90%-decay length for symmetrically and unsymmetrically laminated cylinders, as a 
function of laminate anisotropy. 
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Fig. 3 Effect ofiamina material properties on nondimensional orthotropy parameter for single-ply, homogeneous, 
specially orthotropic laminates (0.2 < v u < 0.35). 



Fig. 4 Effect ofiamina material properties on nondimensional orthotropy parameter for I(±<))) m ],, [(+<t>) 2 m L> an< * 
[(±*j>) m ] T laminates (m = 1, 2, ...). 
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Orthotropy parameter, 


0 


12 A„D tt 
(A n A M - A J „)h* 


1/4 



Fig. 5 Nondimensional orthotropy parameter for typical laminates made of IM7/5260 material. 



Fig. 6 Effect of lamina material properties on nondimensional orthotropy parameter for [(0 2 /±45)J s laminates. 
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Fig. 7 Effect of lamina material properties on nondimensional orthotropy parameter for quasi-isotropic laminates 
and similar unbalanced laminates. 



Fig. 8 Effect of lamina material properties on nondimensional anisotropy parameter for [(+<)>) 2 m l, 
laminates (m = 1, 2, ...). 
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First-order 

correction 

factor, 

e \ 



Fig. 9 Effect of lamina material properties on nondimensional (irst-order correction factor for [(+»!>)*.], 
laminates (m = 1, 2, ...). 


Nondimensional 
90%-decay length, 


d 

/Rh 


£= 0.1 



Fig. 10 Nondimensional 90%-decay length for [(+<t>) 2 m] s laminates made of IM7/5260 graphite-bismaleimide and 
P-100/3502 pitch-epoxy material (h/R = 1/20; m = 1, 2, ...). 
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Fig. 11 Effect of lamina material properties on nondimensional anisotropy parameter for [(0/90/+45 2 ) m ] s 
and [(+45 2 /0/90) m ] s laminates. 



Fig. 12 Effect of lamina material properties on nondimensional anisotropy parameter for [±<j)] T laminates. 
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Anisotropy 
parameter, Q g 
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0.2 
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Fig. 13 Effect of stacking sequence number on nondimensional anisotropy parameter for [(±<|>) m ] T laminates 
made of P-100/3502 pitch-epoxy material. 
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Fig. 14 Effect of lamina material properties on nondimensional first-order correction factor for [±<})] T laminates. 
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correction 
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Fig. 15 Effect of stacking sequence number on nondimensional first-order correction factor for [(±ij>) m ] T laminates 
made of P-100/3502 pitch-epoxy material. 
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Fig. 16 Nondimensional 90%-decay length for [±<(>] x laminates made of IM7/5260 graphite-bismaleimide and 
P-100/3502 pitch-epoxy material (h/R = 1/20). 
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Fig. 17 Effect of lamina material properties on nondimensional orthotropy parameter for [0 p /90 q ] T laminates. 



Fig. 18 Effect of lamina material properties on nondimensional anisotropy parameter for [0 p /90 q ] T laminates. 
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Nondimensional 
90%-decay length, 


d 
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8 = 0.1 



Fig. 19 Nondimensional 90%-decay length for [0 /90 ] T laminates made of IM7/5260 graphite-bismaleimide and 
P-100/3502 pitch-epoxy material (h/R = 1/20). 
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Fig. 20 Effect of lamina material properties on nondimensional orthotropy parameter for [70 p /0 q ] x laminates, 
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Fig. 21 Effect of lamina material properties on nondimensional anisotropy parameter for [70 p /0 ] T laminates. 
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Fig. 22 Nondimensional 90%-decay length for [70 /0 ] T laminates made of IM7/5260 graphite-bismaleimide and 
P-100/3502 pitch-epoxy material (h/R = 1/20). 
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